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Abstract 

A key component of explaining the array of galaxies observed in the Universe is the feedback of active 
galactic nuclei, each powered by a massive black hole’s accretion disc. For accretion to occur, angular 
momentum must be lost by that which is accreted. Electromagnetic radiation must offer some respite in 
this regard, the contribution for which is quantified in this paper, using solely general relativity, under the 
thin-disc regime. Herein, I calculate extremised situations where photons are entirely responsible for energy 
removal in the disc and then extend and relate this to the standard relativistic accretion disc outlined by 
[Novikov fc Thorne! which includes internal angular-momentum transport. While there is potential for the 
contribution of angular-momentum removal from photons to be >1% out to ~10^ Schwarzschild radii if 
the disc is irradiated and maximally liberated of angular momentum through inverse Compton scattering, 
it is more likely of order 10^ Schwarzschild radii if thermal emission from the disc itself is stronger. The 
effect of radiation/scattering is stronger near the horizons of fast-spinning black holes, but, ultimately, other 
mechanisms must drive angular-momentum liberation/transport in accretion discs. 

Keywords: accretion, accretion discs - black hole physics - galaxies: active - galaxies: nuclei - relativistic 
processes - quasars: general 


1 INTRODUCTION 

It has long been widely accepted that active galac¬ 
tic nuclei are powered by gravitationally liberated en¬ 
ergy from accretion discs around massive black holes 
( Lynden-Bell||1969 ). While the study of accretion is an 
interesting prospect in itself, the feedback it ensues is 
significantly consequential for the evolution of galaxies 


as well (e.g. Di Matteo, Springel & Hernquist 2005). 


To truly gauge the effect of feedback requires highly de¬ 
tailed cosmological hydrodynamic simulations that self- 
consistently track the growth of black holes and the 
emission from their accretion discs. However, even the 


most state-of-the-art simulations presently (e.g. Vogels- 
berger et al.||2014 Schaye et al.||2015 ) are unable to re¬ 
solve accretion discs, and must use subresolution mod- 


els to describe their physics (e.g. 

Springel, Di Matteo & 

Hemquist|2005 

Booth & Schaye|2009 

), in a similar vein 

to semi-analytic models (e.g. Croton et al. 2006 

Benson 


20121. An analytic understanding of the functioning of 


accretion discs is hence key for this cause. 


* astevens@swin.edu.au 


In the simple picture of a thin accretion disc, parti¬ 
cles quasi-statically shrink on equatorial, circular orbits, 
where pressure forces are assumed negligible, until they 
reach the orbit of lowest energy, after which they are 


assumed to be captured by the black hole (Lynden-Bell 


1969 [Bardeen 1970). In doing so, those particles must 


be liberated of their angular momentum. Either angular 
momentum is lost through the disc out to higher radii, 
or it is emitted vertically and removed from the disc 
entirely. The latter occurs naturally through thermal 
emission of the disc and through scattering of photons 
if the disc is subject to irradiation from an external 
source. It is of interest then to assess the contribution 
of angular-momentum liberation that photons provide 
in order to better understand the process of accretion 
itself. This paper aims to calculate exactly this, using 
purely general relativistic arguments. 

In Sectionof this paper, relevant mathematical for¬ 
mulae for studying accretion discs are outlined. General 
relativistic calculations are performed based on these 
formulae in Section]^ where limiting cases for the lib¬ 
eration of angular momentum via photons, as well as 
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from the standard Novikov & Thorne (19731 disc, are 


considered. Concluding remarks are provided in Section 

El 


2 MATHEMATICAL FORMALISMS AND 
BACKGROUND 

The unique metric for spacetime around a (non- 
charged) rotating source mass (e.g. a black hole) was 
first discovered by Kerr (1963), usually now written in 


Boyer-Linquist coordinates (Boyer & Lindquist 1967), 


for which the invariant interval is 


- c^dr^ = c^dt^ + ^dr2 + p d0^- 


,2,2, r^ra'^ sin^(9) \ . ^ 


+ a^ + 


sin"( 0 ) dcj)^ 


2rsrar sim( 0 ) 


cdt d(j) \ (la) 


p = r^-\-a^ cos^(0) , A = — r^r + 


(lb) 


where c is the speed of light, = 2GM/c^ is the 
Schwarzschild radius, and a = J/Me is the spin param¬ 
eter (specific angular momentum) of the source of mass 
M and angular momentum J. The metric is stationary, 
axisymmetric, and, as Carter (1968) showed (but see 


also Misner, Thorne fc Wheeler 1973 §33.5), exhibits 
four constants of motionQ Two of these constants are 
the azimuthal and time components of covariant four- 
momentum. Taking the limit r —>■ oo, one finds these to 
be relativistic analogues of energy and azimuthal angu¬ 
lar momentum, usually referred to as the energy and 
angular momentum “at infinity”. Often the “at infin¬ 
ity” is dropped for brevity, and the usual symbols for 
these quantities are used, i.e. 


E = -pt 


L. = 


P<p ■ 


( 2 ) 


When discussing the emission or transport of energy or 
angular momentum in accretion discs (or Kerr geometry 
in general), these are the quantities that are meant. 

Throughout the rest of this paper, most quantities 
will be expressed in a dimensionless form, represented 
by a bar placed on the quantity of interest. For quan¬ 
tities with dimensions of distance, this means normal¬ 
ising to half the Schwarzschild radius, e.g. f = 2r/rs, 
a = 2a/Vg, in line with literature convention. Equation 
[^covers quantities with other dimensions. 

By analysing equations of motion for particles in a 


Kerr spacetime, Bardeen, Press & Teukolsky (1972) 


obtained expressions for the energy and specific an¬ 
gular momentum for circular (i.e. p” = 0 ), equatorial 


(i.e. 9 — 7 r /2 and p® = 0), Keplerian (i.e. gravity is en¬ 
tirely centrifugally balanced) orbits: 


E = 


E 


p3/2 _ 2 ^.i/2 ± a 


TTIC 7-3/4 ^7*3/2 — ^ 2a) 


_xl/2 ’ 


L. = 


L, 


±r^ — 2ar^/^ ± a? 


Vg — 3r^/2 ± 2 a) 


1/2 ’ 


(3a) 


(3b) 


where upper signs are for prograde orbits and lower 
signs retrograde. Note that m represents the rest mass 
of a particle (at infinity), not its inertial mass. By 
checking the derivatives of these quantities. 


dE 8 ar^/^ — 3a^ -I- r(r — 6 ) 
dr W 




(4a) 


- 3 


1-3 - 1 

-a — ar \ -r — 1 


W = 2f^/^ 


2a - 




(r-3) 


1 3/2 


; (4b) 


(4c) 


one finds these two relations share a common minimum, 
referred to in the literature as the innermost stable cir¬ 
cular o rbit (ISCO), where f iscp is given by equation 
2.21 of Bardeen et al. (19721 ^ 

Under the picture where particles transit between in¬ 
finitesimally adjacent orbits in the process of accretion 
(until they reach the ISCO), the above equations pro¬ 
vide the starting point for calculating how much specific 
angular momentum can be lost from photon emission 
and/or scattering. Hereafter, the use of E and (with 
or without bars, but without further subscripts) refers 
to the orbiting states for which these equations 
apply. 


3 RADIATING AND SCATTERING 
AWAY ANGULAR MOMENTUM 

In each of the following subsections, the relative specific- 
angular-momentum loss to photons in an accretion disc, 
—dLz,j/dLz (where subscript 7 is for photons), as a 
function of radius, is calculated for a different idealised 
situation, where each builds on the last. Each result is 
plotted in Fig.j^for r > nscOi allowing for comparisons 
between each individual case. The models considered 
here are all of thin, relativistic accretion discs, whereby 
the mathematics of Section is applicable. 


^In fact, those authors showed t his for the more gener al, charge- 
inclusive Kerr-Newman metric ([Newman et al.||1965|. 
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3.1 Pure, relatively isotropic emission 

Before considering consequences for other forms of en¬ 
ergy transport, the energy carried away by a photon 
supplied by a particle moving to an adjacent lower- 
energy circular orbit should, at most, be the energy 
difference between the orbits. This can be written in 
terms of differentials as 


= —AE . 


( 5 ) 


If this energy is lost primarily through radiation, one ex¬ 
pects each photon to be emitted with statistical isotropy 
from the frame of the emitting particle. For a non¬ 
rotating disc, this would make the average direction of 
emission from each face of the disc vertical (i.e. initially 
completely in the ±0-direction). For a rotating disc per¬ 
ceived from an external frame (i.e. one static with the 
Boyer-Lindquist coordinates), this can then be modelled 
by stating that photons are emitted in the <j)-9 “plane” 
with a three-velocity in the (/(-direction equivalent to 
that of the disc, naturally a function of radius. This (an¬ 
gular) velocity is found as = p^/p^. Recognising = 
—g**E + g^’^Lz a.nAp'^ = —g^'l’E -|- (cf. Equations 

and [^, obtaining the contravariant metric compo¬ 
nents by taking the matrix inverse of Equation and 
expanding and simplifying with Equat ion [3} one con¬ 
cludes consistently with Bardeen et al. (1972) that 


pH- 


±{2rsY^‘^c 


pt 2r^/‘^ ± (2rs)^/^a 


( 6 ) 


Let us write the (contravariant) four-momentum 
components of an emitted photon as Ap^. One can then 
simultaneously solve 


5ttdp* -k gt^Ap'^ = -AE^ 


( 7 ) 


and Equation (for the latter, the left-hand side now 
reads dp'^/dp‘) to obtain explicit functions for dp* and 
Ap'^. Further calculating Ap^ = ALz^-^^ one obtains 


dL 


2,7 


AE-y 


E 


±r^ — 2ar^/'^ ± of 
2 (f3/2 _ 2ri/2 ± a) ’ 


( 8 ) 


consistent with the report of |Lynden-Bell|| ( [i986[ §3.10). 
Now through Equation combined with use of Equa¬ 
tion 1^ one finds an explicit form of —ALz,'^/ALz, as 
presented by the dot-dashed lines in Fig. [T] For non¬ 
spinning black holes, the analytic relation simplifies to 


dL 


^^(a = 0) = 


dL 


f- 2 ’ 


( 9 ) 


which further reduces to the Newtonian case presented 


by Johnson (20111 for r 3> rj. 


Under this picture, one finds that specific-angular- 
momentum removal by photons is important beyond 
the percent level out to ^SOr^. For non-spinning holes, 
as particles approach the ISCO, the radiative efficiency 



Figure 1. (Specific) angular momentum lost to photons relative 
to the (specific-)angular-momentum gap between adjacent, equa¬ 
torial, circular orbits in a relativistic accretion disc as a function 
of radius. Thick curves apply for an accretion disc around a non¬ 
spinning black hole, while thin curves are for a hole spinning with 
a = 0.998 (the maximum of|Th orne |1974| | . The dot-dashed curves 
assume photons are emitted with energy equal to the difference 
of the orbits and are angled to the accretion disc plane such that 
the ^-velocity of the photons matches that of the disc itself (see 
Section |o_ The solid curves follow the solution of [Novikov &:| 
[Th orne |l |1973| l, which inclucie effects of internal torques (outlined 
further in Section |3.2| l. The dashed curves also account for inter¬ 
nal torques, and show the upper limit of scattering, where the 
momentum imparted on photons is parallel to the (/(-direction 
(Section|3.3|l. 


of angular momentum approaches 25%, while for max¬ 
imally spinning black holes, the efficiency approaches 
87%. Already from this analysis, it is clear, and per¬ 
haps unsurprising, that radiation is insufficient by itself 
to liberate an accretion disc of its necessary specific an¬ 
gular momentum. 


3.2 Relatively isotropic emission with 

internal angular-momentum transport 

The previous subsection considered a limiting role 
of photons without any additional form of angular- 
momentum transport. Because photons are unable to 
remove all the necessary angular momentum if they re¬ 
move all the necessary energy, some other mechanism 
must remove angular momentum, which consequently 
must alter the energy liberated by photons too. 

The standard model of thin accretion discs proposed 


by Shakura h Sunyaev (1973), and extended to be rela- 
tivistic by [Novikov Lc Thorne (19731, considers angular 
momentum to be transported radially through internal 
torques (generated by magnetically induced “viscosity” 
- see|Lynden-Beli||1969| [Shakura fc Sunyaeq|1973[ |Lin,| 


Liu h Xiaoqing|20^ ), in addition to removal from pho¬ 

tons, in an accretion disc whose structure is completely 


PASA (2015) 

doi: 10.1017/pas. 2015.XXX 














































4 


A. R. H. Stevens 


stable (i.e. not a function of time). If dLz^nt/df repre¬ 
sents the radial net removal of specific angular momen¬ 
tum from internal torques during a particle transition 
to an infinitesimally adjacent orbit, then it must hold 
that 


dLz ^ dLz,int dLz_^ _ ^ 

df df df 


( 10 ) 


With some small rearranging, the solution derived 
purely from the continuity equations of rest mass, angu¬ 
lar momentum, and energy by Novikov & Thorne (1973) 
provides 


dUzjint 

dr 


(11a) 

dZz.7 

3Zz S 

(lib) 

df 

2 f2 ’ 


where S, S®, and are dimensionless quantities that 
approach unity for increasing r: the reader is referred 
to equation 5.4.1 of [Novikov fc Thorne ( 1973| ) for their 
formal definitions^] 

As for the previous subsection, —dL^ ^y/dL^ is pre¬ 
sented for the Novikov & Thorne (1973) model in Fig.[^ 
as solid lines. For non-spinning black holes, the libera¬ 
tion of angular momentum from photons at low radii is 
noticeably less effective than the previous limiting case. 
Faster spinning holes reach a peak efficiency of radiative 
angular-momentum removal of nearly 40%. 

Naively, one may have expected the solid lines in 
Fig. 0 to lie underneath the dot-dashed lines (i.e. the 
result of Section 3.1), because a new mode of angular- 
momentum transport has been introduced since Section 


3.1 However, as noted by Shakura & Sunyaev (1973), 


internal angular-momentum transport provides a net 
energy source for r ^ risco: meaning more angular mo¬ 
mentum must be liberated by photons, specifically by 
a factor of 3. This is reconcilable by considering an en¬ 
ergy outflow balance equation (i.e. an energy version of 
Equation 10), 


dA dAint _ Q 

dr dr df 


( 12 ) 


Using a combination of Equations Hji, § and 
one can determine —dE^/dE to show that indeed this 
asymptotes to a value of 3. This is plotted in Fig. 
(solid curves). When interpreting Fig.one should ap¬ 
preciate that the absolute energy gap between adjacent 
orbits tends to zero as f —)■ oo, and that thin discs have 
higher densities toward their centres. 


^Page & Thorne| ||l974ll note a sign error for equation 5.4.Ih. 
Those authors also provide an analytically integrated expression 
for SI , but using that provides different results and does not sat- 
is fy Equation|10| wher eas numerically integrating the equations 
of [Novikov &: Thorne] ||1973| | does. For an alternate parametri- 
sation, see [Rhlern^iero 



Figure 2. Relative energy scattered away or emitted through 
photons between infinitesimally adjacent orbits for a |Novikov| 
I&; Thorne] | |1973[ | accretion disc (solid curves, Section El' and 
a maximally scattering-dominant disc with internal angular- 
momentum transport (dashed curves, Section |3.3[ l around non¬ 
spinning (thick curves) and maximally spinning (thin curves) 
black holes. Where the curves pass above a value of 1 (dot-dashed 
line), extra energy is radiated away, transferred internally within 
the disc. The dotted line indicates an asymptote; at large radii, a 
massive particle in the disc emits a photon with thrice the neces¬ 
sary energy to reach its adjacent lower orbit to account for energy 
supplied by internal transport. 


3.3 Scattering with internal 

angular-momentum transport 

If an accretion disc is irradiated by an external source, 
incoming photons can be absorbed or scattered, causing 
particles in the disc to lose angular momentum, a la the 


Poynting-Robertson effect (Poynting 1904 Robertson 


1937 Burns, Lamy & Soter 1979). In the absorption 


case, where the motive absorbing particles re-emit the 
radiation, the analysis of Section |3.2| remains sound. 
However, it could also transpire that charged particles 
in the disc anisotropically transfer their energy to the 
incoming photons via inverse Compton scattering. 

An external irradiative source is observationally mo¬ 
tivated by X-ray reflection spectra generated by an ac¬ 
cretion disc’s corona, which provide a means for measur¬ 
ing black holes’ spins (for a review, see |Reynolds|2014 ). 
At very high redshift, the cosmic background radiation 


could also provide an irradiative source (e.g. 

Fukue & 

Umemura|1994 

Mineshige, Tsuribe & Umemura|1998 ), 


although far more modest in temperature. For fast¬ 
spinning holes, a notable portion of emitted radiation 
from accretion discs is expected to fall back on the discs 


as well (Cunningham 1976). So long as the Thomson 


regime is applicable in the scattering particle’s refer¬ 
ence frame, the scattered photon can have its energy 
significantly multiplied and be beamed in the direction 
of the scatterer’s motion as perceived by an external 
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observer (see Rybicki & Lightman 1979, §7.1). While 
the precise direction of the photon’s change of momen¬ 
tum would depend on its initial energy relative to the 
scattering particle, one can consider the extreme upper 
limit where this is the ())-direction, as shown immedi¬ 
ately below. 

Let us now write the four-momentum components 
of that imparted on the scattered photon as dp^ [it is 
very important to note that this is the change in pho¬ 
ton’s momentum from scattering; when considering ra¬ 
diation, this was the momentum of the (average) photon 
itself, as it did not exist prior (i.e. it had no initial mo¬ 
mentum)]. Because photons have no rest mass, it must 
hold that p^j.dp^dp'" = 0 , where the terms are the 
metric components obtainable from Equation hence 

5 tt(dp*)^ + 25t0dp*dp'^ -f g^^idp'^f = 0 . (13) 

Simultaneously solving Equations and recognising 
= dp 0 = gt4,dp* + g^p^dp'^ in this maximal limit 
(cf. Equation]^, and using Equation]^ one obtains 


dL 


^,7 


dE^ 




±r\/a^ -|- — 2 r — 2 a 


(14) 


One can consider taking the same energy lost to pho¬ 


tons from Section 3.2 but instead using it to kick (scat¬ 


ter) photons in the ())-direction. The relative increase in 
specihc-angular-momentum loss can then be found by 
taking the ratio of Equation [T4| to Equation]^ However, 
by increasing dL^_.y/dr, it must be true that dL^^p^t/dr 
decreases, in order to satisfy conservation of angular 
momentum (Equation [l0|). If dL^p-at/df decreases, then 
so must dEint/dr by an amount found by taking the 
ratio of Equation to Equation after making the 
internal-torque terms the arguments for each: 

dl 

z,int _ dLz/dr+ (dLz,.^/dr)^ 


dE; 


int 


dE/dr + (dE^/drtal ’ 

where subscri pt |3.2| implies the quantities as determined 
from Section Consequently dE^/dr must increase 
from energy conservation (Equation 12), and therefore 
dL^_.y/dr must be higher than initially calculated. One 
can iteratively work through these calculations until 
finding a converged result]^ 

Using the above method, one can calculate the max¬ 
imum —dLz^j/dLz for photon scattering with the ef¬ 
fects of internal transport included. This is shown by 
the dashed lines in Fig. Consistent with the above 
results, under this maximal regime, accretion discs can 
remain efficient above the percent level for —dLz^'y/dLz 
beyond lO'^rg. 

As was the case for the standard INovikov fc Thornel 


(1973) disc, —dE^/dE asymptotically approaches a 


value of 3 for increasing radii, but does so more slowly. 


^One can update the terms with subscript |3.2| in the iterations, 
but the end result is the same. 


It also exceeds a value of 1 at the same radius, as dis¬ 
played by the dashed curves in Fig. 

It should again be stressed that the calculations in 
this subsection are an upper limit. In truth, one should 
expect the relevant curve on Fig. to lie between the 
solid and dashed ones presented, with a bias towards 
the former. 


4 CONCLUSION 

As material accretes onto massive black holes, there 
must be a process by which specific angular momen¬ 
tum is removed from the system. As accretion discs 
are known to be bright sources of radiation, the emis¬ 
sion of photons provides one channel for this angular- 
momentum loss. If a disc is irradiated, inverse Comp¬ 
ton scattering provides another channel. This paper has 
provided calculations of the contribution of angular- 
momentum liberation through photons in thin, rela¬ 
tivistic accretion discs. In addtion to situations where 
transport by internal torques was included, such that 
the usual conservation laws of physics were satisfied, 
these calculations included limiting situations where 
photon emission was responsible for all the necessary 
energy removal and where photon scattering was an¬ 
gled to remove maximal angular momentum. 

On scales up to the order of a hundred Schwarzschild 
radii, photons remove a small (> 1 ) percentage of an¬ 
gular momentum in accretion discs that emit in an ex¬ 
pected fashion, but beyond this contribute negligibly. 
At the absolute most, discs subject to strong irradiation 
are potentially capable of scattering away angular mo¬ 
mentum as efficiently out two orders of magntiude far¬ 
ther from the black hole. The contribution in both cases 
becomes stronger near the horizon of a fast rotating 
black hole, especially in the latter case (<60%, cf. thin 
dashed line. Fig. [^. By and large though, angular mo¬ 
mentum is transported far more efficiently through the 
disc internally, rather than liberated from the disc elec- 
tromagnetically. 
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